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Abstract. We propose a model of magnetic-field growth in galaxy clusters whereby the field is amplified by a factor of about 
10* over a cosmologically short time of ~ 10** yr. Our model is based on the idea that the viscosity of the intracluster medium 
during the field-amplification epoch is determined not by particle collisions but by plasma microinstabilities: these give rise to 
small-scale fluctuations, which scatter particles, increasing their effective collision rate and, therefore, the effective Reynolds 
number. This gives rise to a bootstrap effect as the growth of the field triggers the instabilities which increase the Reynolds 
number which, in turn, accelerates the growth of the field. The growth is explosive and the result is that the observed field 
strength is reached over a fraction of the cluster lifetime independent of the exact strength of the seed field (which only needs 
to be above ~ 10~^^ G to trigger the explosive growth). 
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1. Introduction 

Rapidly improving measurements of Faraday Rotation in 
galaxy clusters reveal that intracluster medium is permeated 
by tangled magnetic fields with strength of a few and 
typical coherence scales of 1 — 10 kp c ( iKr onberg. 1994; 

Li 



iCarilli & Taylor , 



Clarke. Kronberg & BohringeJ. l200 it . , , 

20Q2t iGovoni & FerettiL bOOllVogt & EnBlinI l2Q03i 12005). 

This means that the intracluster medium (ICM) contains 
huge amounts of magnetic energy and it is natural to ask 
whether there is a robust and universal mechanism that 
could be responsible for generating such fields and that 
would also be insensitive to the particular circumstances in 
individual clusters (specific age of the cluster, magnitude of 
the seed field, etc.). We would like to forego the discussion 
of whether cluster fields could be purely of primordial or 
external origin (as siiggested by several authors; see, e.g., 
IKronberg etall l200lHBaneriee & Jedamziil '2003') and in- 
stead show that they can be amplified very quickly and with 
no stringent constraints on the magnitude of the seed field by 
the small-scale dynamo effect due to the turbulent motions 
of the ICM. The idea that the fields a re dynamo generated 
has been explored by many authors jjaffel Il98(i iRolandl 



1981 


iRuzmaikin. Sokoloff & Shukiiro^ 


19891 iDe Young 


1992 


iGoldshmidt & Rechaelil Il993l 


iKulsrud et al. 


1997 


, Sanchez-Salcedo. Brandenburg & Shukurov' Il998 


.Subramanian. Shukurov & Haugen. 2005 


: EnBlin&Vogt 


I2OO5 


). The main issues are whether there is, indeed, turbu- 



lence in clusters and whether the rms rate of strain associated 
with this turbulence is sufficiently large to exponentially 
grow the field to the observed strength in the cluster lifetime 
of a few Gyr. 

There are many possible energy sources in clusters that 
almost certainly produce random fluidlike motions of the 
ICM; merger events, subcluster and galactic wakes, active 
galactic nuclei. Vario us indire ct observa tional estimate s (e.g. , 
ISchueckeret a l.L Eo"o4; .Rebusco et al.[ 12005- Fuiitl Eooi 
appear to converge in expecting random flows with veloc- 
ity fluctuations U ^ 10^ — 10'^ km/s at the outer scale 
L ^ 10^ — lO'^ kpc (a direct detection may be achieved 
in the near future; see llnogamov & Sunvaevl [20031) . Similar 
number s are obtained in n umerical simulations of cluster for- 
matio i i jNorman & Brvar .>1999:.Sunvaev. Norman & Brvanl 
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l2003t iRicker & Saraziil l200lHTakizawaL l2005b . There is. 
however, no consensus on whether turbulence, at least in the 
usual hydrodynamic sense of a Kolmogorov energy cascade 
across a broad inertial range of scales, is a generic feature 



1 INTRODUCTION 



of clusters The main difficulty is the 

very large values of the ICM viscosity obtained via the stan- 
dard estimate ~ wthiAmfp, where vthi ~ 10'^ km/s is the 
ion thermal speed and Amfp ^1 — 10 kpc is the mean free 
path (assuming fully ionised hydrogen ICM with tempera- 
ture T ^ 10® K and density n ^ 10^^ cm~'^). This gives 
Re ~ 10^ if not less, which makes the existence of a well- 
developed inertial range doubtful. 

In fact, the small-scale dynamo does not require a turbu- 
lent velocity field in the sense of a broad inertial range. As the 
dynamo effect is simply the random stretching of field lines, 
it is always controlled by the rms rate of strain — and this is 
associated with the viscous-scale motions, which are random 
but spatially smooth. Whether the viscous scale (~ LRe"^/"^) 
is much smaller than the outer scale L is inessetial as long as 
random motions are maintained' — not a problem in clusters, 
given the variety of available stirring mechanisms. Thus, the 
field will be amplified and the only issue is whether this hap- 
pens sufficiently fast, given the relatively small growth rate 
obtainable in a viscous ICM. We shall see that this turns out 
not to be a problem because the effective ICM viscosity in 
the weak-field regime is much smaller than WthjAmfp. 

The widely used viscosity estimate ~ wthiAmfp only ap- 
plies when the magnetic field is so weak that that the ion gy- 
roradius is larger than the mean free path: pi > A^fp. This 
gives B < 10^'® G — much smaller than the observed field 
in clusters. For stronger fields, the plasma becomes magne- 
tised and its viscosity anisotropic: the standard value still 
works parallel to the magnetic field, while the motions per- 
pendicular to the field are virtually undamped (Brasinskii, 
[l96 5). As was pointed out recently by ISchekochihin et alJ 
( l2005i) . the anisotropy of the viscous stress in clusters means 
anisotropy of the plasma pressure, which, in turn, always 
gives rise to well-known plasma microinstabilities (most im- 
portantly, mirror and firehose, see, e.g.. iGarvi.il993ii) at scales 
all the way down to the ion gyroscale. The emergence of the 
pressure anisotropics is intimately related to the changes in 
the strength of the magnetic field: in a magnetised plasma, the 
conservation of the first adiabatic invariant for each particle, 
fj, = v\/B, implies that the perpendicular plasma pressue 
changes according to^ 

1 dp_L I AB pj_- p\\ 



(1) 



p± dt H dt p± 

where the first term is due to the conservation of p, and 
the second term represents collisional isotropisation (uu ^ 
t^thi/Amfp is the ion collision frequency). Balancing the two 
terms, we get 



,;A, 



1 dB 

where A = {p± — p^^)/p±. On the other hand, since the re- 
sistivity of the ICM is tiny (Rm ^ 10^^), the field Unes are 

' See numerical simulations of 
with Ro ^ 1-10^ by ^ 
iHau een. Bra ndenburg & DobleJ 



(2) 



small-scale dynamo 
Schekochihin et al.' (■2004'); 



^004) and the review of the 

astrop hysical turbulence and dynamo by ISchekochihin & CowlevI 
<2005l) . 

^ Here and in the rest of this paper, really means 



frozen into the fluid flow and the change of the magnetic-field 
strength can be expressed in terms of the fluid velocity u: 



1 dB 
B~dt 



= bb 



Vu-V ■u= ^Re'/^ 



(3) 



where b = B/B. The last equality in Eq. Q is based on 
estimating the turbulent rate of strain at the viscous scale. It 
is the parallel viscosity that matters here because the motions 
that change B are precisely the o nes — and the only ones — 
that are damped by it ( Bra ginskiil 1 1 965h . Since 



Re = 



UL 



UL 



we have, from Eqs. (|3 and Q, 

1/2 / \ 2 



A = 



U 



U 1 



U 

Vthi 



Re 



-1/2 



(4) 



(5) 



The peak growth rate of these instabilities (at the ion gy- 
roscale) is 



2 V 

7,nax = ( |A| - - I Vti 



(6) 



where (3 = SirnT/B^, fl^ = eB / crrii is the ion cyclotron fre- 
quency, ai — \ for the mirror and ai — 1/2 for the firehose 
instability. To sum up briefly, external energy sources drive 
random motions, which change the field [Eq. (|3|l], which 
gives rise to pressure anisotropics [Eq. (|2j], which trigger the 
instabilities. The latter are stabilised if the magnetic field is 
sufficiently strong: (3 < 2/ A. 

It is not currently clear exactly how the instabilities affect 
the structure of the turbulence and magnetic fields in clusters. 
One plausible hypothesis is that the saturation is quasilinear 
with small fluctuations SB at the gyroscale: 



SB^ 
B^ 



2X "2 



(7) 



where a2 is some positive power. These fluctuations will scat- 
ter particles, giving an effective collision frequency: 



^scatter 



SB^ 



■In 



(8) 



where a = ai + 0:2 > 0. The result is an effective Reynolds 
number [Eq. @ with vu I'scattor] that is much larger than 
the one based on the actual ion-ion collisions and that in- 
creases with the magnetic-field strength (because I'scattcr c>c 
rii oc B). Since this, in turn, speeds up the small-scale dy- 
namo by increasing the turbulent rate of strain [Eq. Q], we 
expect explosive growth of magnetic field. This po ssibility 
was briefly referred to in ISchekochihin et al .' (2005). In §|2| 
we shall construct a simple model of such a self-accelerating 
dynamo and evaluate its effectiveness in generating observed 
cluster fields for a set of fiducial cluster parameters (intro- 
duced in §|2ji. We will conclude with a summary and some 
words of caution in §|3 



3 A MODEL OF SELF-ACCELERATING CL USTER DYNAMO 



3.1 Exponential stage 



2. Cluster parameters 

In order to makes specific estimates as we proceed, let us 
assume a fiducial set of parameters for the ICM:^ 

n = IQ-^ cm^^ (9) 

T = 10** K ^ vtM = \ — ^ 1000 km/s, (10) 

V 

U = 300 km/s (rms turbulent velocity), (11) 

L = 200 kpc (outer scale), (12) 

Amfp = 1 kpc. (13) 

We shall measure magnetic field in units of 

Seq = (87rnT)i/2 ~ 20 fj.G, (14) 

for which P{Bcq) — 1. For this field, the ion cyclotron fre- 
quency and the ion gyroradius are 

^i,cq — — 0.2 S , 



crrii 

Vthi 



a 



5000 km. 



(15) 
(16) 



The magnetic field strength above which plasma is magne- 
tised (pi < Amfp) is 



Bo — Be 



A 



10-^» G. 



mfp 



(17) 



The turbulence is characterised by Re ~ 60 and 



h ^ yr 
L 



(outer timescale), (18) 



U 



Re 



-1/2 _ inS 



10 yr (viscous timescale), (19) 



~ LRe ^ IQ kpc (viscous cutoff). 



(20) 



3. A model of self-accelerating cluster dynamo 

We assume that the net effect of the instabilities is to change 
the effective collision frequency of ions: 

^cS = Vii + Vsci,ttcAB), (21) 

where i^scattor is given by Eq. (|8|i. If we use Eqs. ( 1415 1 with 
Vu replaced by v^s, we find the following equation for the 
anisotropy ratio 



where magnetic field is in units of i?oq [Eq. ( I14H and 

1/2 

- 0.01, 



Ao 



U 



A 



mfp 



u 



Pi.cq _ Pi.cq « 2 



L 



A 



10 



-17 



mfp 



(22) 



(23) 



(24) 



Our results will not be very sensitive to the value of a. In 
what follows, we shall use a — 3/2 for specific estimates. 

' We emphatically do not indend to suggest that the actual values 
of these quantities are known for real clusters to more than an order- 
of-magnitude precision (in some cases, even worse than that). These 
numbers are merely a convenient example. 



As the instabilities grow much faster than the field 
changes, we can assume that, as far as the field is concerned, 
the adjustament of the collision frequency occurs instanta- 
neously. The growth of the magnetic field satisfies [from 
Eqs. Q and Q with i^u i^cs] 

-- = ^ (25) 
B dt A{B) ' ^ ' 

where time is in units of [vthi/UY'L/U and A(_B) is the 

solution of Eq. M2\ . Equations (I22> and M5\ constitute our 

model of the field growth. They are solved anlytically in three 

asymptotic regimes in what follows and numerically in Fig.fD 

3.1. Exponential stage 

Let us start with an arbitrarily small seed field. As long as 
the the first term in Eq. (I22> dominates, the effect of the in- 
stabilities is negligible and the field grows exponentially with 
the characteristic e-folding time ^ Aq, which in dimensional 
units translates to the inverse rms rate of strain ^ 10* yr 
[Eq. ( I19H . This initial dynamo stage continues until the sec- 
ond term in Eq. M2\ becomes comparable to the first, which 
occurs at the field strength 



Bi = 



. 2+Q 



= B, 



Pi, cq 



cq 



mfp 



For a = 3/2, we have Bi 



U J Amfp 

10-15 G. 



a/2 



(26) 



3.2. Explosive stage 

For B ^ Bi,we may approximate the solution of Eq. ( I22t by 

/ £ X l/(2+a) 

MB) ^ (I) . (27) 

Upon substitution into Eq. (I25t . this gives an explosively 
growing field: 

where, taking -8(0) = Bi, we get tc = {2 + a)Ao. During 
this very fast stage, the anisotropy drops precipitously, while 
the effective Reynolds number and, therefore, the rms rate of 
strain, grow: 

Rcoff - ( —] -T^- (29) 



B' 



,l/(5+2a) 



B, 



(30) 



^VtMj A(B)2- 

Hence the self-accelerating nature of the dynamo. Eventually 
(in fact, very soon), the field grows to be strong enough so 
that B^ is comparble to A{B), the second term in Eq. i22\ 
starts dropping again and the A stops decreasing and starts 
increasing. The minimum value of A is reached at B ^ B2, 
where 

, l/(5+2a) 

U \ Pi^ 
Vthi J L 

For a = 3/2, we have B2 ^ 10^^ G. The time it takes to 
achieve this is 

1 ~ - (2 + a) e2/(5+2") ^tc^ 10* yr. (31) 
Over this time, the field is amplified by a factor of ~ 10*. 
This is to be compared with the conventi onal dynamo (e.g., 
ISubramanian. Shukurov & Haueenl l2005l) . which has the e- 
folding time of 10* yr and would, therefore, need around 

2 Gyr to achieve an amplification factor of 10*. 



3.4 Nonlinear stage 
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(a) (b) 




Fig. 1. (a) The anisotropy ratio vs. field strength. Dotted lines are the analytical solutions for the three asymptotic regimes 
discussed in §|3lll §ll2l and §E3l[A = Aq, A(S) given by Eq. and A{B) = 2B^]. (b) Growth of the magnetic field. 
The dotted lines show various field-strength thresholds discussed in the text [Eqs. (I17> . i26\ . (I30> . ( I34> and ( I36H . The initial 
slope of the curve is the rate at which the conventional dynamo can grow the field. Both plots are for a — 3/2 and the fiducial 
set of cluster parameters given in §|2l 



3.3. Algebraic stage 

When B > B2, A stays just above 2B^. Computingthe small 
correction, we obtain the following solution of Eq. i22\ 



A{B) ~ 2B^ 



1/0 



B 

Solving Eq. (|25} with A ~ 25^ and 5(0) = B2, we get 



(32) 



B{t) = \Bl+t 



(33) 



This is a rather slow growth but it only lasts until the first 
term in Eq. \22\ is again dominant: the transition is, in fact, 
sharp because once A < 2B^, the instabilities are shut down 
and we have to set i^scattcr = for all greater values of B. 
The field at which this happens is 

-,1/4 



Ap 
2 



1/2 



B, 



eg 
V2 



U 



mfp 



1 nG. (34) 



The time to get there is t ~ i?| = Ao/2 ~ 10^ yr again. 
The amplification factor since the end of the explosive stage 
is - 10. 

Note that that if we compare pi with the effective mean 
free path A,„fp,off = Vt\,i/vcS, we get 

P^{B) ^ -Beg £ 
A„fp,cff(S) B A{BY' 

which is readily seen to decrease throughout all of the dy- 
namo stages discussed above. Therefore, the magnetised- 
plasma assumption, once true at _B = Bq, is never broken. 
For example, at the point of the largest effective Reynolds 
number (Re ~ 10® at end of the explosive regime, B = B2) 
we have Amfp,off — 0.03 pc and pi ~ 10® km. 

3.4. Nonlinear stage 

So far, the dynamo we have considered has been kinematic, 
i.e., the magnetic field has remained too weak to exert a 
back reaction on the flow. The condition for this to be true 



is that the field energy remains smaller than the energy 
of the smallest-scale turbulent motions that can change the 
field, i.e., the viscous-scale motions (see*Scheko chihin et all 
|2D04). In Kolmo gorov turbulence, their velocity is [/ Re ', 
whence, using Eqs. ( I29l l and ilAl . we find that the condition 
for the nonlinearity to become important is B^ > A{B). 
Thus, all through the algebraic regime, the magnetic field 
hovers around this threshold. After instabilities are stabilised, 
it can finally exceed it by growing above B^. For B > B3, we 
have to model the effect of the back reaction. How to do this 
is a highly nontrivial problem for which a definitive solution 
has not so far been found. Here, we shall use the model pro- 
posed by Schekochihin et al. (2002),"^ which is constructed 
in a somewhat similar spirit as the self-accelerating dynamo 
model above. The main idea is that the growing magnetic 
field gradually suppresses the ability of turbulent eddies to 
stretch it. At any given point in time, all eddies whose energy 
is less than the energy of the field are thus suppressed. Fi- 
nal saturated state is reached when the energy of the field is 
comparable to the energy of the outer-scale eddies, i.e., to the 
total energy of the fluid turbulence: 



B4 = Seq — ~ 6 fiG 

Vthi 



(36) 



To model this process, we assume that the magnetic field 
evolves according to 



1 dB 
B~dt 



(37) 



with the S-dependent growth rate dSchekochihin et alll2002l) 

1 - 



U 



7(5) = ^Rc^/^ 



1 



[l + Bl/Blf 



1 



{l + B^/BlY {I + Bl/BIY 



1/2 



(38) 



For numerical evidence in support of the model, see 
.Schekochihin et aLt.2004i) . 
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It is not hard to see that, for i?3 < i? < Bi, -f{B) cx 
which again gives the magnetic field growing as 
Thus, the slow algebraic growth continues as the field ex- 
ceeds -B3 until it saturates at B = B4. Again, it does not take 
long because the inertial range in clusters is very narrow and 
the field does have to grow very much to get from B^ to B4. 

4. Conclusions 

All stages of the evolution of the magnetic field described in 
the previous section are illustrated in Fig.^ which presents 
the numerical solution of Eqs. i25\ and i22\ [followed by 
Eqs. ( I37I38> in the nonlinear stage]. We see that the field 
grows by a factor of about 10^^ over a time period of ^ 
2 Gyr. Most of the growth occurs over a short fraction of 
this time during the explosive self-accelerating dynamo stage, 
which takes the field from any value above ~ 10^^^ G to 
^ 10^^ G in a cosmologically short time of about 10^ yr. 
After that, the field does not have to grow further for a very 
long time before it reaches the strength of several /iG consis- 
tent with the observational data. The importance of the self- 
accelerating dynamo mechanism we have proposed is that 
it provides for field amplification to the observed strength 
over a short period of time independent of the precise value 
of the seed field^ or, indeed, of exactly how many Gyr a 
cluster has been around. This removes the need to envi- 
sion ways of generating seed fields satisfying the often quite 
large lower bounds that theories (cited in the Introduction ) 
and numeri cal simulations (e.g.. | R pettinger. Stone & Burns. 
199^ J)olag. Bartelmann & Le s'dTTl999.. .2002; .Dolag et al.. 
2005c iBriiggen et al.. 2005) infer to be necessary for mag- 



netic fields of observed strength to be produced by the 
conventional turbulent dynamo over the cluster lifetime. 
Since the amplification time is so short, it is also proba- 
bly not crucial whether the cluster turbulence responsible 
for amplification is best modelled as forced or decaying (cf. 
ISubramanian. Shukurov & Haugen, 2005). We conclude that, 
no matter what the seed field is, the random motions in the 
ICM will have no difficulty in amplifying it to the observed 
level in a fraction of the cluster lifetime. 

Finally, we must alert the reader to an important is- 
sue, which we have ignored above. This concerns the 
spatial structure of the cluster fields. Theory and simu- 
lations of the small-scale dynamo in the MHD descrip- 
tion indicate that the dynamo-generated fields have a 
folded structure with direction reversals at the resistive 



scale ~ rail99d 


Schekochihin et all l2004 


iBrandenburg & Subramanianl 


I2OO5I 


). In clusters, this scale 



turns out to be extremely small (~ 10^ km), while analysis of 
the observa tional data suggests that the t ypical reversal scale 
is ~ 1 kpc jVo gt & Enfiliil l2003[ Eool . It is likely that the 
reversal scale is, in fact, set by some form of anomalous re- 
sistivity associated with electron scattering by the plasma in- 
stabilities discussed above. We relegate further development 
of this idea to future work. 



For seed fields below ~ 10^ G, the conventional exponen- 
tially growing dynamo with e-folding time ~ 10® yr is required 
until 10^^^ G is reached. 
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